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Abstract 

We present a direct approach to the construction of Lagrangians 
for a large class of one-dimensional dynamical systems with a simple 
dependence (monomial or polynomial) on the velocity. We rederive 
and generalize some recent results and find Lagrangian formulations 
which seem to be new. Some of the considered systems (e.g., motions 
with the friction proportional to the velocity and to the square of the 
velocity) admit infinite families of different Lagrangian formulations. 
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1 Introduction 

In recent papers [H [2] a problem of finding Lagrangian description for a 
large class of one-dimensional dissipative (or dissiptive-looking) systems was 
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discusssed. The discussion was far from being exhaustive. In this paper we 
present a different, more direct, approach to the problem of the construction 
of Lagrangians for disipative (or dissipative-looking) systems. We simply 
assume some general form of the Lagrangian and then check the resulting 
Euler-Lagrange equations. 

The inverse problem of Lagrangian mechanics is concerned with the ques- 
tion whether a given system of second-order ordinary differential equations 
= q-, q) can be derived from a variational principle [3j. In other words, 
one tries to find a Lagrangian for this system. This problem was studied in 
XIX century by Helmholtz (see [1]) and by Darboux who proved that in the 
one-dimensional case the Lagrangian always exists [S]. The inverse problem 
in two-dimensional case was solved by Douglas [B], while the general case has 
been completed recently 

Dissipative systems were long believed to be 'beyond variational treat- 
ment' |8], which is to some extent true if we insist on physical interpretation 
of the Hamiltonian and canonical momenta, compare [9]. However, by re- 
laxing these requirements one can get variational interpretation of numerous 
dissipative system pUl [TTl [T^ [T3] . 

In this paper we focus on more elementary issues, namely on providing 
explicit Lagrangian description for a large class of one- dimensional differential 
equations of second order with a simple (e.g., polynomial) dependence on the 
velocity x. 



2 Standard Lagrangians 

Standard Lagragians (known also as 'natural' or 'of mechanical type') are 
quadratic forms with respect to x (the dot denotes the differentiation with 
respect to t). In the one-dimensional case we can easily obtain all equations 
of motions corresponding to standard Lagrangians. We assume 

£ = ip(x,t)i;^ + g(x,t)x + i?(x,t) , (1) 
The Euler-Lagrange equations yield 

i + §i^ + ^i + ^^=0, (2) 

where subscripts x, t denote partial derivatives. As a consequence we imme- 
diately obtain the following proposition. 
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Proposition 2.1 The equation of motion 

X + a{x, t)x^ + t)x + c(x, t) = (3) 
admits a Lagrangian description with a standard Lagrangian (QP iff 

K = 2at . (4) 
Then P = exp(2 a{^,t)d^) , and 

R = j\Qt{^,t)-ci^,t)Pi^,t))d^ , (5) 

where Q = Q{x,t) is an arbitrary function. 

All systems described by ([1]) have also the Hamiltonian description. In- 
deed, computing generalized momentum 

p = Px + Q , X = ^ (6) 
we easily get the standard Hamiltonian H = px — C: 

i^(a:,P,t) = ^^^^^-i?(x,t). (7) 

Corollary 2.2 Special cases of Proposition [2A[ 
1. P = P{t) andQ = 0: 

This is a generalization of Proposition 1 from JE/- In the case of linear 
equations (i.e., c = xc{t)) we have: 

In particular, we rederive the well known result [T4\ for the damped 
harmonic oscillator: 

X + jx + uJqX = =^ C = -e'^^ (i^ — cUqX^) . 
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2. P = P{x) and R = 0; 



This formula simplifies for c = c{x) (the case considered in JUl^). 
Then 

C= f^i;2 + ti;c(x)V'^'"^^^'^ (9) 



^2 

3. P = P{x) and Q = 0; 

x+a{x)x^+c{x,t) =0 =^ C = ^x^e^^^''^^^^^-J\{^,t)e^^'''^y^''yd^ 

This is a generalization of the main result of ^1] and Proposition 3 from 
l^, where c = c{x). Thus these results are extended on t-dependent 
function c = c{x, t). 

4. P = A{x)B{t): 

X + a{x)x'^ + b{t)x + c{x,t) =0 =^ L is given by ([I]), where 

P = AB, A = exp(2 ra{^)d^), B = exp(/* 6(r)rfr), R is given by 
and Q is arbitrary. 



Example 2.3 (A particle accreting mass in a potential field) We pro- 
ceed to physical aspects of the equation 

x + b{t)x + c{x,t) = . (10) 

Following ITB^ . where the damped harmonic oscillator is interpreted as har- 
monic oscillator with time- dependent mass, we define 



(t)=e/*^W^^ i.e., &W = ^- (11) 
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Then, 

C = im(t)i2 _ mit)Vix, t) , H = + m(t)\/(x, t) , (12) 

where V{x,t) = c{^,t)dC,. Therefore, the equation [IW can be consid- 
ered either as a dissipative system or a particle with a prescribed mass time- 
dependence in an arbitrary potential (possibly time dependent) . 



3 Reciprocal Lagrangians 

Reciprocal Lagrangians (i.e., inverses of standard-like Lagrangians) were in- 
troduced and studied recently [21 [171 HH] • If 



£ = Y, L = L{x,x,t)), (13) 



then 



2^dLdL ^j^d'^L ^OLdL ^d'^L ^dL 
■■ _ dx dx dxdx dt dx dtdx dx 

dx"^ \ dx 

We confine ourselves to L of the form 

L = F(x,t)i'^ + G(x,t) . (15) 

Substituting f[T^ into ffT^ we obtain: 

px^^ + qx^"'^ + rx^ + sx^~^ + w , . 

gx''-^ - hx^'^-^ ' ^ ' 

where 

p ■.= {! + v)FR,, q := vFFt, r := (1 + 2v)FG^ + (1 - v)F,G, 

s := 2i^GtF - uGFt, w:=GG^, g:=v{v-l)FG, h:=v{v+l)F^. 

The case z/ = 1, F = 1 is discussed in [17J, with a special stress on G 
quadratic in x (leading to second-order Riccati equations), see also [2]. In 
the case z/ = 1 the equation flTB]) reduces to a special case of ([3]): 

.2 (Fi + 3G,). {2GtF-GFt + GG,) 

x = — —X X — . 17) 

F 2F 2F2 ^ ^ 



First, we confine ourselves to t- independent F, G. Then the coefficients 
a, b, c by powers of x depends on x only. They are not independent. Indeed, 



F' 3G' GO' 

a = — , = , c = , (18) 

F 2F 2F^ V y 

where the prime denotes the differentiation with respect to x. Hence, sub- 
stituting G = 3cF/b and F' = aF to the last equation of f|T8|) . we get a 
constraint on a,b,c, see f|T9|) . 



Proposition 3.1 The equation ^ admits the Lagrangian description with 
L = (if (i) + G(i))-i iB 

Then, F{x) = exp(J^ a{^)d^) and G{x) = 3c{x)F{x) /b{x) . 

Therefore, we can choose arbitrary functions a{x),b{x) and then c have 
to satisfy the equation (IT^ . Solving this equation we get: 

c{x) = h{x) J 6(0 exp (^J^ a{y)dy^ d^ . (20) 



Example 3.2 Taking a = and b{x) = kx (k = constj we obtain 

c{x) = \x (hzx^ + = ^ + Aix , (21) 

where A = const and \\ := ^kX. This case corresponds exactly to a Lienard- 
type nonlinear oscillator which shows very unusual properties, like isochronous 
oscillations for Ai > 118f . 

Another possibility is to choose arbitrary functions b{x),c{x) and then a{x) 
is given by: 

a = ^ - ^ + — . (22) 

b c 9c ^ ^ 



Proposition 13.11 generalizes Propositions 4 and 5 from [2]. The case u = 1 
contains other interesting subcases. Indeed, assuming F = f(t), G = xg(t) 
the equation f|T71) can be reduced to the hnear equation: 



X + h{t)x + c{t)x = , (23) 

where 

2f ' 2/2 • ^ ^ 

The system expresses b, c in terms of /, g. It turns out that these equa- 
tions can be inverted. Given b, c we may compute corresponding /, g. Indeed, 
substituting g = 2/6 — / into the second equation we get the inhomogeneous 
linear equation: 

which can be solved in quadratures in a standard way: 

fit) = exp (^1* (26(r) - c(r)) (^exp £ b{y)dy^ dr^ dz^ . (26) 



Proposition 3.3 The equation [2B\) (for any b{t), c{t) ) admits a Lagrangian 
description with the reciprocal Lagrangian of the form C = {xf{t) + xg{t))~^ , 
where f is given by and g = 2/6 — /. 



Therefore, any equation of the form fl23|) (including equations of mathe- 
matical physics, like Airy, Bessell, Hermite or Legendre equation) admits at 
least two different Lagrangians: standard one (see (jHD) and reciprocal. 



Example 3.4 We get another simple case taking F{t) = foC and G{t) 
goc^^. Then [Tl\) reduces to x + kx = 0. 
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In the case v = 2 the Lagrangian fllSI) yields more comphcated equation: 

.. _ 3FF.x" + 2FFtx' + {hFG, - F,G)x^ + {AGtF - 2GFt)x + GG^ 

^ 2F(G'-3Fi;2) ^^^^ 

In the particular case G = G{t), F = f{x)G^ we obtain: 

^ = l^i' + . (28) 



Proposition 3.5 The equation 

X + a{x)x'^ + h{t)x = (29) 

admits a Lagrangian description with a Lagrangian given by [Fx^ + G)^^ , 
where G{t) = exp(- /* 6(r)dr), F{x,t) = exp(-3 f b{T)dT - 2 fai^dO- 

Therefore, the equation fl29l) admits at least two different Lagrangian 
descriptions: standard (compare the case 4 of Corollary 12. 2p and reciprocal. 

4 Lagrangians with a modified kinetic term 

In this section we consider generalizations of standard Lagrangians, where 
the kinetic term x^ is replaced by some more general expression (and the 
term linear in x is absent). 

4.1 Monomial case 

First, we assume the monomial case: 

C = F{x,t)x^' -G{x,t). (30) 
The equation of motion reads 

^~ 'ilF' ~ (/X - l)F ~ /i(/i - l)F ■ ^ ' 
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Proposition 4.1 The equation of motion 

X + a{x,t)x'^ + b{x,t)x + c{x,t)x'^'^ = (/i 7^ 0, 1) (32) 
admits a Lagrangian description with the Lagrangian ^3^) iff 

{fi - l)b,cc = fia,t . (33) 
Then F = exp(/i f a{^, t)d^) and G = - 1) c(^, t)F{^, t)d^. 

The proof follows directly by comparing (1321) with (1311) . Another result 
is obtained by assuming F = F{x) and G = G{x). 

Proposition 4.2 The equation 

X = -a{x)x'^ - c(x)x^ [v ^ 1, 2) (34) 
admits for any a(x),c(x) a Lagrangian description. The Lagrangian reads 

C = F{x)x^~'' - G{x) , (35) 

where 

Fix) = exp (^(2 - u) j\{Od^^ , G(x) = (2 - z/)(l - u) J\{OFim- 



Corollary 4.3 Taking c{x) = 0, a{x) = k = const, and denoting n = v 
we obtain (for n ^ 0) 

x + kx^ = =^ C = Ci;"e"*^^ . 
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4.2 General case 

Let us consider a class of standard-like Lagrangians with quadratic kinetic 
terms replaced by an arbitrary smooth function of x. 

C = t)i;{x) + G{x, t) . (36) 

The equation of motion reads 

(Ft + xFM'-F^tfj-G^ 
^ + — 'l^. ^ = 0. (37) 

Assuming F = Fq = const we get the equation 

where the right hand side is of the form f{x,t)(j){x) for some functions /, 0. 
Indeed, it is enough to take Gx = fFo oraz ip" = 1/0. 



Proposition 4.4 The equation x = f{x,t)R{x) admits a Lagrangain de- 
scription with the Lagrangian L = "^{x) + G{x,t), where 

^{v) := 1' dr^l'^^, G{x, t) = 1^ /(e, t)d^ , (39) 

(provided that the above integrals exist). 



Corollary 4.5 Special cases of Proposition 4-4 



1. X = xf{x,t) =^ £ = xln|i;| + J f{^,t)d^ 



3. x = -kox'' =^ C = - --kox (ly^ 1,2) 

(2 - z/)(l - u) 

4. x = f{x,t)(l-^] =^ c = ~c\l-%+ f{tt)di 



(?■ j V 



10 



5 Radical Lagrangians 

We consider Lagrangians of the form 



L = ^A{x,t)x'' + B{x,t). (40) 



The Euler-Lagrange equations yield 

_|_ ^j,u _ 
gx2i^-2 _|_ 



px'^" + qx'^'^ ^ + rx'^ + sx" ^ + w 



where 



(l-/i) 
z/(l-^) 



(U - Ufl- fl) + 1) 




fi B^B 

' * (l-/i) ^ ' ' z/(l-/i) A 



(42) 



In this paper we will assume either /x = z/7^1or/i7^z/ = l. In those 
cases the denomiator simplifies and the righ-hand side of (14TI) is a polynomial 
in X. 

5.1 The case ^ = u ^ 1 

In this case the equation ( l4Ti) reduces to 



a; 



2Ai^;^+2 _|_ At -iz+l _|_ / (l+y)A^ _ uB^\ ™2 , / {l-i')Bt , iMt \ • , B^-2~u 
"r^-^ "r| A B I ^ \ B ~^ A I ^ A 



v{l - V) 

(43) 

A further reduction is obtained by assuming that A = A{t), B = B(t). Then 
the equation becomes 



Proposition 5.1 The equation 

X = -a{t)x - b{t)x''+^ (45) 

admits (for u ^ 0, u ^ 1 and any functions a,b) a Lagrangian description 
with the Lagrangian of the form C = \J A{t)x^ + Bit) where 

rt - - Nl-- 



A{t) = (ur h{r) exp {-u f a{y)dy) dr 

(46) 

(-Z/ a{y)dy) dr^ exp (^-u f a{r)dT^ . 

In order to proof this proposition it is enough to compare (H5ll with (j4 
and to solve resulting differential equations. 

Two interesting special cases can be obtained by requiring either 6 = 
(i.e., A{t) = const) or a = (i.e., ulnA — {u — l)lnB = const). 

Corollary 5.2 Special cases of Proposition [3771 - 
1. A = Ao = const, B = B{t): 



X + a(t)x = 



C= a Aox" + Boexp{-u j b{T)dT) 



where Bq = const and u ^ 0,1. 
2. B = CqA^, 

X + h{t)x'^ = =^ C= + cqF-^-i 

where m 1,2, F = F{t) = — co(m + 1) /* h{T)dT, and cq = const. 
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5.2 The case 7^ z/ = 1 

In this case the equation fHTl) reduces to 

.. _ + fi)x^ + {At + i?.. (1 - 2/x) + x + Bt{l-fi) + i^-i^ 

(47) 

We assume A = A{t), B = B{t). Then the equation ( 1471) becomes 
■■ - ^ ■ ME _ B 

and, solving linear differential equations (similarly as in the case of Proposi- 
tion [5?T]), we get the following result. 

Proposition 5.3 The equation 

X = a{t)x + b{t) (49) 

admits (for any functions a,b) a Lagrangian description with the Lagrangian 
of the form C = A{t)x + B{t) where /i 7^ 1 and 

A{t) = exp ^(/i — 1) J* a{T)dT^ , 

(50) 

B{t) = - (^/*6(r)e-i'""(s')'^^rfr) exp (^fi f a{T)dT^ . 



We point out that for 6 = formulas (|50l) yield B(t) = exp ( /i f a{r)dT 
(the integration constant has to be taken into account). 

6 Multi-Lagrangian cases 

The Lagrangian of Proposition 15.31 can be rewritten as 



jr = e/* <^(^)d^ ^ xe- /* <^)<^^ - j h{T)e- f ^y'^'^y , (51) 

and this form suggests the following generalization which can be easily veri- 
fied by a simple straightforward calculation. 
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Proposition 6.1 The equation 

X = a{t)x + b{t) (52) 

admits (for any functions a,b) a Lagrangian description with the Lagrangian 
of the form 

..e/-.-^(..-r"--/'.Me-r--). (53, 

where F is a function of one variable (such that F" 0). 

In particular, we point out tliat tlie simple classical equation x + kx = 
has Lagrangians of all forms considered in our paper, namely: 

Li = -e'^'x^ , L2 = ^^-^ TT , L3 = x^e^^'-^'>''' , 

2 e^^^x + e'^* ^g^-j 



L4 = i; In — kx , L5 = ^/ x'^ + e 



ukt 



However, most of these forms (except L4) can be reduced to the Lagrangian 
of Proposition 16.11 namely 

Lf = e-^'F{xe''' + cq) , (55) 

with F{S,) equal to v^, respectively. 

Another multi-Lagrangian case is described by Corollary 14.31 where we 
present a one-parameter family of Lagrangians for the equation x + kx^ = 0. 
What is more, the corresponding Hamiltonian is proportional to the La- 
grangian (for any n ^ 0) and is time-independent. Hence, C is an integral of 
motion. This observation can be generalized as follows. 

Proposition 6.2 Suppose that a Lagrangian C = C{q^,q^,t) is an invariant 
of motion (i.e., dC/dt = 0). Then, for any (sufficiently smooth) function 
F : M. ^ R, the Lagrangian C = F[C) yields the same equations of motion. 

The proof is straightforward. We compute: 

dC, _ dF dC dC _ dF dC d dC _ d'^F dC dF d dC 
dqi~dZdq^' d4'~dCW' '(UW' ~ d£^^ ^ dCdiW ' 
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Therefore, 

d dC _dC _ f d dC _dC\dF d'^F dC 
diW ~W~ [diW ~W)dC^ d£^~dt ' 

from which the proof follows immediately. 

Taking into account Proposition 16.21 we see that Cp := F{xe^*) is a La- 
grangian for the equation x + kx^ = (for any smooth function F) . An- 
other Lagrangian (time- independent) for this equation was found by Sarlet: 
C = x{l — Inx) exp{kx), see |19j|. 

7 Conclusions 

In this paper we succeded to rederive all results of PQ [2] in a straightfor- 
ward, simple way. Actually, we found many other dissipative-looking sys- 
tems possessing a Lagrangian description. One- dimensional systems admit- 
ting the Lagrangian formulation were discussed in numerous papers (see, e.g., 
[l6l[I7l[l8l|19l|20l|2ll|22l|23llM[25l|2^ some of them devoted mostly to 
the damped harmonic oscillator, e.g., [TH [T51 EZl- Surprisingly enough, us- 
ing quite elementary tools, we succeeded to find some number of Lagrangians 
which seem to be overlooked in the exisiting literature. 

It is interesting, that for some systems the Lagrangian description is not 
unique: they may belong to several classes (the problem of the equivalence 
was first discussed in [28]). The equations x + kx = and x + kx"^ = 
0, usually considered as classical dissipative equations (compare [S]), have 
infinite families of Lagrangians, see Section [61 The equation x + kx = has 
Lagrangians of all forms considered in our paper, see fl5^ . 
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